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A broad range of quantum optimisation problems can be phrased as the question whether 
a specific system has a ground state at zero energy, i.e. whether its Hamiltonian is frustration 
free. Frustration-free Hamiltonians, in turn, play a central role for constructing and under¬ 
standing new phases of matter in quantum many-body physics. Unfortunately, determining 
whether this is the case is known to be a complexity-theoretically intractable problem. This 
makes it highly desirable to search for efficient heuristics and algorithms in order to, at least, 
partially answer this question. Here we prove a general criterion - a sufficient condition - 
under which a local Hamiltonian is guaranteed to be frustration free by lifting Shearer’s 
theorem from classical probability theory to the quantum world. Remarkably, evaluating 
this condition proceeds via a fully classical analysis of a hard-core lattice gas at negative 
fugacity on the Hamiltonian’s interaction graph which, as a statistical mechanics problem, 
is of interest in its own right. We concretely apply this criterion to local Hamiltonians on 
various regular lattices, while bringing to bear the tools of spin glass physics which permit 
us to obtain new bounds on the SAT/UNSAT transition in random quantum satisfiability. 
These also lead us to natural conjectures for when such bounds will be tight, as well as to 
a novel notion of universality for these computer science problems. Besides providing con¬ 
crete algorithms leading to detailed and quantitative insights, this underscores the power of 
marrying classical statistical mechanics with quantum computation and complexity theory. 


An overwhelming majority of systems of physical interest can be described via local Hamilto¬ 
nians: 


M 

^ = En. (1) 

i=l 

Here, the ‘/c-local’ operator Hj acts on a A:-tuple of the microscopic degrees of freedom, best thought 
of as qudits for the computer scientists among our readers, or spins for the physicists. The M 
operators define an interaction (hyper-)graph G, as displayed in Fig. 

A surprisingly diverse and important class of such model Hamiltonians is defined by the addi¬ 
tional property of being frustration-free', the ground state \'if) ol H \s a simultaneous ground state of 
each and every Hj. This class comprises both commuting Hamiltonians - for which [Hj, Hj] = 0 Vi, j 
- such as the toric code, general quantum error correcting codes and Levin-Wen models [IHS]; as 
well as non-commuting ones, such as the AKLT and Rokhsar-Kivelson modelsjlHG]. Their par¬ 
ticular usefulness is also related to the fact that many of these examples can be viewed as ‘local 
parent Hamiltonians’ for generalized matrix product states [7]. In general, frustration-free condi¬ 
tions provide analytic control of ground state properties in otherwise largely inaccessible quantum 
problems. 

Determining whether a given Hamiltonian H is frustration-free is well known in quantum com¬ 
plexity theory as the quantum satisfiability (QSAT) problem. QSAT is widely believed to be 
intractable, in the sense that no general purpose classical or quantum algorithm can efficiently 
determine whether a given Hamiltonian is frustration-free (‘satisfiable’). The technical statement 
is that QSAT is QMAi-complete [H], even when restricted to qubits and /c = 3 [S], or when the 
interactions is between neighboring qudits on a line cni. 

Fortunately, such hardness results only apply in the worst case. For instance, in contrast to 
the hardness result, it is immediately obvious that a fully disconnected interaction graph G can be 
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analysed efficiently by considering each term Ilj individually. 



Figure 1. Left: The interaction graph of a fc = 3-local Hamiltonian. The degrees of freedom, qudits, are 
denoted by circles; the squares, which indicate the operators H^ appearing in H, are joined to the qudits on 
which they act. In the hard-core lattice gas mapping, each square may be occupied by a particle covering 
also the adjacent circles (shaded blobs), on which it must not overlap with another particle. Right: A fc = 6 
interaction graph forming a triangular lattice of the operators H^. The corresponding hard-core lattice gas 
corresponds to the hard hexagon model, which is exactly solvable. 

The central result reported here is a sufficient combinatorial criterion for a local Hamiltonian 
H to be frustration free. In fact, we provide a lower bound for the dimension of the satisfying 
subspace. This amounts to a generalisation of Shearer’s theorem|ll] from classical probability 
theory to the quantum world. 

We first formulate the result in Theorem 1, followed by an intuitive explanation of its content, 
with a technical proof relegated to the supplemental information. We then turn to applying 
Theorem 1, for which we enlist the results available on statistical mechanics of hardcore objects with 
negative fugacity on the interaction graph to deduce statements regarding QSAT, producing new 
bounds on the satisfiability threshold for a large class of one, two and three dimensional interaction 
graphs. With help of the cavity method, we are able to conjecture improved lower bounds for the 
satisfiability of QSAT on regular and Erdos-Renyi random graphs, canonical models for quantum 
constraint optimization problems. These statement hold just as well for classical satisfiability - 
but for some of these classical models better bounds are known nans]. 

We close with an outlook, including a discussion of the role of a universality which emerges in 
our analysis, as well as conjectures on when our results are exact or tight. 

Theorem 1. Given a local projector Hamiltonian H as in Eq. 0 with interaction graph G and 
relative projector rank p = R(Imnj) for all i, then 

R(kerH) > Z{G,-p) > 0 

if Z{G, —p') > 0 for all 0 < p' < p. 

Without loss of generality, we have taken the terms Hj to be projectors so that the satisfiability 
condition reduces to the presence of a zero energy ground state |14j . 

Here, the relative dimension of a subspace X of the full Hilbert space H is simply R(W) = 
and, Z(G, A) is the partition function for a hardcore lattice gas of fat particles living on the (hyper)- 
edges of the interaction graph at fugacity A (see Fig. [^. More precisely, 

Z{G,X) = Y.X^^^l[{l-n,n,) ( 2 ) 

{m} i<r^j 

where i ^ j runs over all projectors which share qudits and the sum runs over occupations Ui = 0,1 
of the lattice gas. 
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Intuition — Let us first consider the case of classical projectors, Ilj, diagonal in the standard 
computational basis, where we can provide a simple pictorial representation of the content of the 
Shearer bound. The operator HiCl “blj) projects onto the zero-energy space of the Hamiltonian. If 
it is non-zero the Hamiltonian is frustration-free. Moreover, its expectation value in the maximally 
mixed (infinite temperature) state gives the relative dimension of the zero-energy space, 

/2(kerH) = n(l-ni) (3) 

If all of the projectors act on different qudits, then, this expectation value simply factors, 

M 

i?(ker H) = n (1-n.) = (1 - p)^ 
i=\ 

where p = Hj is the relative dimension of Hj. This is non-zero so long as p < pc = 1. 

Now suppose that some of the projectors share qudits; we call such pairs of projectors dependent, 
as they introduce dependence into the expectation value in Eq. (§. When can we, nevertheless, 
guarantee that R(keT H) remains positive? The essence of the Shearer bound is captured by the 
Venn diagram in Fig. where each bubble represents the fraction of configuration space ‘knocked 
out’ by a projector, with overlapping bubbles representing configurations multiply penalised by the 
corresponding projectors. 





Figure 2. Shearer bound for dependent projectors. Projectors are independent if they do not share any 
qudit (top left), and there will always exist configurations violating any combination of such projectors 
simultaneously. In the Venn diagram (top right), this is represented by the mutual intersections of the 
shaded circles, each of which denotes the fraction of configuration space penalised by the projector it is 
labeled with. By contrast, for dependent projectors (bottom left), such an overlap is not guaranteed. A 
lower bound for the relative dimension of the satisfying space can be obtained by assuming the projectors 
Hi and H; do not share any violating configurations with the projector 11^ with which they share a qudit. 
(bottom right) 


The area covered by the bubbles is largest if they do not overlap, corresponding to a large 
number of violating configurations; we thus expect a lower bound on R(kei H) when they are fully 
disjoint. To calculate this lower bound, we expand the product in Eq. ^ and introduce a collection 
of occupation variables n* = 0,1 which indicate the presence of Hj in each term of the expansion, 

i?(ker H) = ^ (-1)^ (4) 

{rii} 

If two dependent projectors, Hj and Hj, are thus occupied, then that term is zero when their 
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bubbles are made disjoint. Otherwise, it is given by Thus, 

R(kei H) > 11 (1 - rurij) 

{rii} i-^j 

= Z{G,-p) (5) 

where i ^ j runs over projectors that share qudits. 

We have derived inequality ([^ under the assumption that it is possible to make the dependent 
bubbles disjoint. If p is small enough, this is always the case; in the Venn diagram, the bubbles 
can be made disjoint without covering more area than the total space contains. Shearer miiEi 
showed that the above intuitive lower bound is correct for classical projectors so long as p < pc 
where Pc is the first zero of Z(G, —p). This is the classical analog of our quantum generalization. 
Theorem 1. 

The classical sketch above makes little sense for non-commuting quantum projectors. In the 
language of probability, this reflects the failure of the inclusion-exclusion principle (Eq. @) for 
the relative dimension of vector spaces. Nonetheless, the result holds; the proof - our fundamental 
technical advance - is provided in Appendix |X} 

Statistical mechanical transcription — A remarkable aspect of Theorem 1 is that it maps the 
quantum satisfiability problem onto the classical statistical mechanical problem of determining 
the position of the first negative fugacity zero \c{G) of the partition function Z{G, A). In general, 
evaluating Z is computationally hard|16]. In the thermodynamic limit (number of qudits V —)• oo), 
the first zero Xc{G) can be identified with a well known critical point Xc{Goo) of the hard-core lattice 
gas referred to as the hard-core singularity. The critical fugacity Ac(Goo) upper bounds the Ac of 
any finite subgraph by monotonicity m- Thus, Theorem 1 implies that all subgraphs of Goo are 
frustration free so long as the relative rank of the projectors satisfies p < — Ac(Goo)- 

The hard-core singularity has been studied extensively in the statistical mechanical literature 
[TTHSS! and its location is known for many infinite lattices. Table |l] summarizes some of these 
critical values and indicates their translation into QSAT interaction graphs. In addition, there 
exist rigorous bounds on Ac which may be efficiently calculated on any sufficiently simple lattice 
m- The singularity exhibits universal features just like standard phase transitions. In particular. 


Table I. Summary of the critical threshold for various infinite interaction graphs 


Lattice 

Pc 

Qudits 

Projectors k 

1-D chain 

i m 

vertices edges 

2 

Triangular 

UnjimiMj edges 

vertices 

6 

Square 

0.1193 [IH1I23I 

edges 

vertices 

4 

Square 

0.0889 im 

vertices edges 

2 

Checkerboard 0.0688 jUj 

black 

red 

4 

Hexagonal 

0.1547 im 

edges 

vertices 

3 

Simple Cubic 

0.0744 US] 

edges 

vertices 

6 

f-regular tree 

mm 

vertices edges 

2 


the free energy density /(A) = ^ logZ{X) near the critical point Ac has universal exponents [211123] 
due to its connection with the so-called Yang-Lee edge singularity HailSlESHSI]. Technically, this 
means that the leading non-analytic part of the free energy / ~ (A — Xc)^^ near Ac where (j)D 
is a non-integer critical exponent which depends only on the spatial dimension D of the lattice 
[22] l23l [25] . It is known analytically that (pi = 1/2, p 2 = 5/6 [HI EH ES] and po = 3/2 for D 
above the upper critical dimension Dc = 6 [25] . 
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The critical exponents control the lower bounds on the dimension of the ground space of 
Hamiltonians close to the critical threshold. Indeed, by Theorem 1, R(kei H) > oc 

gan(A-Ac)'^o+o(A-Ac) £qj. ^ gQ critical exponents show up directly - and perhaps somewhat 
unexpectedly - as anomalous growth laws in Z) = 1, 2. In the following, observation of the expected 
exponent in the random Z-QSAT ensemble helps confirm the validity of the non-rigorous cavity 
analysis. 

Examples — As a first application of Theorem 1, we calculate the critical threshold Ac for 
several infinite graphs. We do this using the cavity method, a well-known technique for studying 
statistical [M] and quantum j35H37| models on infinite graphs which are locally tree-like. On trees 
and chains, the results are rigorously exact while on infinite random graphs they are less rigorous, 
but often just as exact. 

The heart of the cavity method is to introduce an auxiliary system of messages which propagate 
in both directions between the hyperedges and the nodes of the interaction graph G. For the model 
given by Eq. (§, the derivation is straightforward, and quite intuitive at positive fugacity where the 
cavity messages correspond to marginal distributions for the gas particles - see Appendix [B| The 
resulting equations continue smoothly to the negative fugacity regime. Adopting the convention 
that indices i, j,... label hyperedges and a, b ,... label sites, the belief propagation (BP) equations |5H] 

Qa—>-i 




where and Qa^i give the probability that hyperedge i is occupied in the absence of the con¬ 
nection to a and in the absence of everything except a respectively. We note that a special case 
of these equations has been used previously to study the number of dimer coverings of random 
graphs {k = 2, positive fugacity) [39] . 

Given a solution of the BP equations, the cavity estimate of the free energy F = logZ of the 
system is given by a sum of the free energies associated to the addition of individual elements of 
the interaction graph, i.e.. 



^ = + ( 6 ) 

® i {ai) 

where Fa, Fi and Fat correspond to the change in free energy due to the addition of sites a, 
hyperedges i and the links ai between them, respectively. The full expressions may be found in 
Appendix 

The ehain — On the infinite chain, the BP equations can be solved by uniform messages 
qa^i = q, h^a = I- After some algebra and taking the root of the BP equations corresponding 
to q{\ = 0) = 0, one finds I = q = 1 + This expression suggests that Ac = —1/4 as g 

becomes complex for A < Ac. Indeed, it is easy to check that the free energy density f = F/N 
Eq. Q has the expansion near Ac: 

/ = — log 2 -|- 2(A — Ac)^^^ — 2(A — Ac) -!-••• 

The free energy goes complex for A < Ac, which reflects the accumulation of partition function 
zeros, and indicates the failure of the lower bound. 
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The free energy density is precisely — log 2 at the critical point. For qudits of local dimension 
q, this provides a lower bound on the dimension of the ground space dim(ker H) = R(kei H) > 
( 2 )^■ In fact, careful application of Theorem [l for finite chains agrees with the exact result 
dim(ker if) = (|)^ {N + 1) for open chains derived using matrix-product techniques [271140| . It is 
interesting to note that the entropy per site is zero for qubits at criticality, but has a finite value 
for larger q. 

Regular trees — The BP equations on infinite regular trees with site degree t and hyperedge 
degree k can be solved under the ansatz that all messages are the same. For algebraic details, see 
Appendix The resulting critical threshold is 


1 (/ c - 1)^-1 

P 


(7) 


For k = 2-local trees, this agrees with previous results obtained by matrix product states on 
trees j28]. As far as we are aware, Eq. Q provides a strictly better lower bound on satisfiability of 
infinite regular trees than any previous literature. The corresponding expansion of the free energy 
density is 


/ = hit, k) + A(t, k){\ - Ae) + B{k k){\ - A,)3/2 + . 


( 8 ) 


We thus recover the expected non-analyticity {(p = 3/2) in the free energy density for the infinite¬ 
dimensional hard core singularity j25] . 

Random k-QSAT — We now apply Theoremto random f-QSAT. Random satisfiability has 
been a workhorse for the study of typical case complexity and heuristic algorithms. By tuning the 
density a = ^ of Erdds-Renyi-type random interaction graphs, a cornucopia of phases and phase- 
transitions in the structure of the satisfying space has been discovered. In the quantum case, three 
phases for random f-QSAT on qubits and rank-1 projectors are believed to be: 1) a low density, 
a product state satisfiable phase (PRODSAT); 2) at intermediate densities for k sufficiently large, 
a satisfiable phase in which all ground states are entangled (ENTSAT); and 3) a high density, an 
unsatisfiable (UNSAT) phase. The arguably most interesting of these phases, ENTSAT, has been 
shown to exist using the Quantum Lovasz Local Lemma (QLLL) [4T] for k > 13. The application 
of Theorem 1 to the Erdds-Renyi (ER) ensemble is not straightforward, as the interaction graphs 
exhibit an unbounded degree distribution, so that a strict application of Shearer’s theorem provides 
no useful information. However, as explained in Appendix such local degree fluctuations do not 
appear to be the source of unsatisfiability in /c-QSAT for qubits. 

Neglecting this local rare-region effect, we calculate the self-consistent solutions to the disorder- 
averaged cavity equations using standard methods (population dynamics) [38]. Eor A > Ac, this 
numerical technique converges to a population of messages which represents the distribution of 
q and I in the infinite random graph, and from which we can directly estimate the disorder- 

df 

averaged free energy / and occupation number density (n) = 7 —-—. The latter is particularly 

d (log A) 

useful, because as A —)■ Ac, we expect (n) to exhibit a square root singularity by universality, i.e. 
(n) ~ (A — Ac)^/^. Observing this singularity allows us to estimate Ac accurately and confirms the 
validity of the numerical approach. By fitting these singularities, we extract a new and improved 
(lower) but non-rigorous threshold for the existence of an ENTSAT phase, k > 7, see Eig. 

Conelusion and open questions — By extending the classical Shearer theorem to quantum me¬ 
chanical systems, we have provided a new statistical mechanical criterion for determining whether 
a local Hamiltonian is frustration free. We have applied this criterion to a large class of regular 
and random Hamiltonians. These instances cover many of the geometries which have been studied 
in quantum complexity and as parent Hamiltonians for wavefunction-based many-body physics. In 
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2 4 6 

Projector density (a) 


Figure 3. Our current understanding of the phase diagram of the random fc-QSAT. For small a, the 
instances are PRODSAT, with the transition out of PRODSAT approaching ap —?> 1 from below as k 
grows [1^. For large a, the instances are guaranteed to be UNSAT, with the best known upper bound|13] 
for the satisfiability transition at Og = 0.573 • 2^ at fc > 4. In between, there may be an ENTSAT phase 
if as > ap. The best previous lower bound for the satisfiability transition, ag was given by the Quantum 
Lovasz Local Lemma (QLLL, up triangles! |41j . and is roughly exponential at large k. This work obtains a 
better lower bound through the Shearer criterion (squares), lowering the threshold for the existence of an 
ENTSAT phase from fc = 13 to fc = 7. 


the context of random satisfiability problems, we have provided a set of new lower bounds on the 
existence of the ENTSAT phase. In particular, these bounds suggest that the ENTSAT phase is 
eminently more reachable in simulations than previously established. 

Theorem 1 depends on the dimension p of the projectors Ilj relative to the Hilbert space dimen¬ 
sion, and not the absolute local dimension q of the qudits. At small q (e.g. qubits), there are many 
Hamiltonians where i?(ker H) is strictly larger than the bound of Theorem 1 (cf. the discussion of 
random /c-QSAT). Nonetheless, for large q we conjecture that Theorem 1 becomes tight. This is 
in sharp contrast to the classical case where the length 4 cycle (periodic chain) already provides a 
counterexample to the analogous statement [H] . However, it is easy to show numerically that this 
counterexample breaks down for quantum projectors. 

If, indeed. Theorem 1 is tight, there are several striking consequences. The geometrization 
theorem |l0] states that R(kerH) is minimized by almost all choices of the Hi’s. Coupled with 
tightness, the lattice gas partition function Z then provides a complete characterization of quantum 
satisfiability for almost all Hamiltonians with large enough qudits. It also directly lifts the univer¬ 
sality of the lattice gas critical exponents to the counting of the ground state entropy of almost 
all quantum Hamiltonians in the frustration free regime. In this sense, the conjecture amounts to 
an even larger scope of transferring insights from classical statistical mechanics into the quantum 
complexity domain. 

While Theorem 1 can guarantee the existence of zero energy states, it does not construct 
them. This is in contrast to classical SAT and commuting QSAT, where efficient constructive 
classical (under the Lovasz and Shearer’s assumptions) and quantum algorithms (under the Lovasz 
assumption) are known |44H47| . 

Analogously constructing the wavefunctions corresponding to the solutions of the non-commuting 
quantum problem would represent a milestone for quantum complexity theory. 
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Appendix A: Proof of Theorem 

Definition 1. For a subspace X <Z V, let R{X) = • We say that X is mutually R- 

independent of the subspaces Yi,..., C F */VS' C [n], R{X Hies Yf) = R{X)R{f},^sy^)- 

Definition 2. The subspaces Xi,..., Xn have R-dependency graph G = ([n], E) if Xi is mutually 
R-independent o/ where r(i) = {j|(i,j) G E}, r'''(z) = r(z) U {z}. 

An independent set (also known as a stable set) of a graph G is a set of non-adjacent vertices. 
Indep(G) is the family of all the independent sets in G. The independent set polynomial is 

I{G,x)= Y. 

SSlndep(G) 

We use the shorthand I{G) = /(G, —p) whenever p is clear from the context. 

Lemma 2 ([l5]). Fix {G,p). The following properties are equivalent: 

1. For every 0 < p' < p, I{G, —p') > 0. 

2. For every indueed subgraph F of G, I{F, —p) > 0. 

If this is the ease, we say that {G,p) has Shearer’s property. 

Theorem 3. Let Xi,... ,Xn subspaees ofV with dependeney graph G, such thatMi G [n], R{Xi) > 
I — p. If {G,p) has Shearer’s property, then ^ I{G,p) > 0. 

For comparison, the following is the classical analogue by Shearer Un- 

Theorem 4. Let Bi,..., Bn be events with dependeney graph G, sueh that'll G [n] Pr(Si) > l—p. 
U [G,p) has Shearer’s property, then Pr(njg[n] Bi) > I{G,p) > 0. 

Proof of Theorem We will first introduce some notations. For S' C [n], we use the shorthand 
^s = n i^g Xi, where Xf/f = V. The subgraph induced hy U C V{G) is denoted G[U]. 

Lemma 5. For any S C T C [n], 

RjXs) ^ R{Xt) 

/(G[S]) - /(G[T]) 

The lemma completes the proof of the theorem, using S = 0 and T = [n]: = 1. 

We prove the lemma by induction on the size of S and assume that T = SU{1}. Let A = S\F(i). 
By partitioning the independent sets in T to ones that contain i and the ones that do not, we have 

I{G[T])=I{G[S])-pI{G[A]). (Al) 

We get a lower bound similar to the previous equation 


R(At) = R(Xr(i) n W n Aa) 

> R(Ar(i) n Aa) + R{Xi n Aa) - R(Aa) 

= R(As) + R(A0 R(Aa) - R(Aa) 

>R(As)-pR(Aa) (A2) 

The first inequality follows from Fact the second equality follows from the fact that A* is 
R-independent of {Xj}j^A, and the last inequality follows from the assumption R(Aj) > 1 — p. 
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Fact 6. 


R{x n y n z) > r{x n z) + r{y n z) - r{z). 


Factfollows from the following property by assigning A = Xr\Z,B = Yf\Z-. 

dim (74 n i?) = dim(A) + dim(i?) — dim(A + B) 


for any two subspaces A, B C V, where A + B = {a + b\a G A,b G B}. 
To complete the proof of the lemma, 


R{Xt) R{Xs) ^ R{Xs)-pR{Xa) R{Xs) 
mr]) ms]) - ms]) -pi{G[A]) ms]) 
_ p{I{G[A])R{Xs) - R(X^)/(G[5])) 
mS]){I{G[S]) - pI{G[A])) 
pI{G[A]) f R(Xs) R(Xj) \ 

I{G[S]) - pIiG[A]) [l{G[S]) I{G[A])) 

^ pI{G[A]) f R(Xs) _ R(Xj} \ ^ 

I{G[T]) \I{G[S]) I{G[A])) - ’ 


where we used Eq. (Al) and (A2| in the first inequality, Eq. (Al) in the third equality, and the 
second definition of Shearer’s property (see Lemma and the induction hypothesis in the last 
inequality. 

□ 


Proof of Theorem [2 This is an application of Theorem]^ Eor a QSAT instance IIi,..., 11^, with 
interaction graph F we define Xi = ker(nj), and G has m vertices, and an edge between i and j if 
Ilj and lij act on a shared qudit. G is indeed an ii-dependency graph for Xi,, Xm- this follows 
from the property dim(A 0 B) = dim(A) dim(R) |41j . 

To complete the proof we observe that I{G, —p) = Z[F, —p) (in graph theoretic terminology, 
this is equivalent to the statement that for any hypergraph F, = Mp^x), where L(F) is 

the line graph of F and Mp{x) is the matching polynomial of F, see [57] ). □ 


Appendix B: The cavity derivations 

The cavity method is a collection of techniques for calculating the properties of statistical 
mechanical models defined on locally tree-like lattices, i.e. on graphs without short loops. Here, 
we apply it to evaluating the free energy F = log Z of the hard-core gas of particles living on the 
hyperedges i of the interaction graph G, as defined by Eq. (|^. 

We introduce a pair of cavity distributions (‘messages’ or ‘beliefs’), pi^a{ni) and pa^i{ni), for 
each link between a hyperedge i and node a. These represent the marginal distribution of the 
occupation variable n* in the absence of its link to site a and in the absence of its link to sites 
other than a, respectively. On a graph with no loops, these distributions are simply the normalized 
partial sums of Z when evaluated from the leaves inward towards the link between a and i. 
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Assuming that the cavity distributions coming into a given node or edge are independent (as 
they are on graphs with no loops), the 2Mk cavity messages satisfy the local relations, 


Pa^i (nj) = 


jeda\i 


. \keda / 


Pi^a {ni) = A"* n 

bedi\a 


Pb^i jni) 

\ni 


A"' n 

jeda\i 


where p' denotes an unnormalized distribution. Properly normalized, each cavity distribution can 
be parametrized by one real number, which we take to be the probability of occupation (recall that 
Ui G {0,1}), 


Qa^i — Pa—>-2(1) li^a — P2^a(l) 

Using this parametrization leads directly to the BP equations in the main text. 

On graphs without loops, the BP equations admit a single global solution for the cavity messages 
qa^i and which corresponds to an exact evaluation of Z. From this solution, one obtains the 
total free energy Eq. ®) as a sum of local contributions due to the free energy of each hyperedge 
Fi, site Fa and link Fai- These free energies are functions of the messages coming into the relevant 
piece of the interaction graph. Explicitly, 


Fa = log 

= log 
Fi = log 
= log 
Fai = log 
= log 


{rij}, j 
jeda 


n (1^ ^ (1 Ik- 

jeda jeda keda\j 

Pb^i (rii) 


E^”‘n 

. bedi 




n (1 - Qb^i )+A 


Qb^i 


.bedi 


bedi 




(1 I'i^a) (1 Qa^i) H“ 


k—^aQa—^i 

A 


This evaluation of the free energy is exact for finite trees. Eor infinite graphs with loops of divergent 
girth, it often produces asymptotically exact results. 

Thus, the main task in evaluating the free energy (and corresponding Z) for a given infinite 
graph is to obtain solutions of the infinite set of BP equations for the cavity messages. 


1. Solving the cavity equations: regular graphs 

Eor infinite t-regular /c-local trees/regular random graphs, the BP equations admit a solution 
with uniform messages Qa^i = q and U^a = P After a small amount of algebra, the self-consistent 
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equations can be reduced to finding the roots of 

x’^ - - X{t - 1) =0 

where 

i = i + {‘-i) (pi^) 

= Ab-‘ - 1) 


The discriminant of Eq. (Bl) is (with z = A(t — 1)) 


(Bl) 


A = (-1)(2) 


1 -1 

0 ••• 

0 

—z 

0 

0 \ 

0 1 

-1 0 


0 

-z 0 ••• 

0 

0 

... 1 

-1 

0 

0 

—z 

k 1 — k 

0 ••• 

0 

0 

0 

0 

0 k 

1 - A: 0 


0 

0 

0 

0 


k 

1 - k 

0 

0 / 








fc+i 



k-2 



= (-1)(2)(-z)^-2 ^-zk^ + (-1)^ ((1 - k)'^ + k{l - 
= {{k - 1)^“^ + k’^z) 

where in the second step we observed that only 3 terms in the determinant are non-zero. The 
discriminant has a simple zero at Zc = — —, in agreement with Eq. ([^. The zeros of the 

discriminant indicate that values of 2 ; where two roots of the BP equations merge so that there is 
a non-analyticity in the resulting free energy. 

At Zc, the root of interest has multiplicity 2 and can thus be easily found by looking for a root 


shared with the derivative of Eq. (Bl |. Thus, 


k-1 


Xc = 


k 


which can be verified by direct evaluation. Geometrically, we expect x to evolve as a square root 
of deviations in z from criticality. To leading non-vanishing order in 6z = z — Zc and 5x = x — Xc, 


Eq. (Bl) can be solved. 


6 z = ^{Sxf 




(B2) 


Erom the critical values for Zc and Xc, and the behavior near criticality Eq. (B2), it is a straight 


forward if algebraically tedious process to evaluate the free energy density given in Eq. 


2. Solving the cavity equations: random graphs 

In infinite Erdos-Renyi random graphs, the local geometry fluctuates from site to site and, 
accordingly, the local BP equations vary. We do not expect to find a solution with uniform 
messages; rather, the distribution of cavity messages Pq[-] and Pi[-] ought to be stable under the 
iteration of the BP equations. This leads to self-consistent equations for the distributions Pq[-] and 
Pi[-], which we solve by numerical population dynamics. 























14 


Appendix C: Local degree fluctuations for random A:-QSAT: 

In an Erdos-Renyi random graph, the degree z of each site follows a Poisson distribution, i.e., 

Pr (degree = z) = e~^" . 

This means that there is a small but finite density of large degree sites at any a. These large degree 
sites pose a problem for the Shearer bound: the partition function of the finite ‘star’ of hyperedges 
surrounding a site of degree z is Z = 1 + z\. This has its first negative fugacity zero at Ac = —Ijz. 
Given the unbounded fluctuations in z, monotonicity implies that Ac = 0 for the infinite graph. 

However, we believe that it is not these fluctuations which in fact lead to the instance becoming 
UNSAT for the case of qubits. This can already be seen classically: one can explicitly verify that 
a star which by Shearer should be UNSAT will, in fact, always be SAT; this happens because 
the constraints imposed by different projectors on the central spin, can never be disjoint as the 
worst-case Shearer bounds assumes. 

This carries over to the quantum case. Indeed, here it has been shown [l3] that stars are not 
only satisfiable, but even have entropic ground spaces dimkeri^ = 2(2^“^ — -|- 1). This 

is backed up by the rigorous work of m, who find a lower bound on the satisfiability threshold 
by gluing together product states on high degree regions of the Erdds-Renyi graph with states 
guaranteed by the QLLL on the low degree regions. 

Our numerics, in turn, also implicitly neglect the large degree fluctuations because the density 
of sites with z > 2^ is very small (10“^'^ at A: = 5 to at A: = 9) for a ~ of the observed 

bulk transition. By contrast, for k < 5, the large degree fluctuations are not as rare and thus the 
population dynamics fails to converge. 

Eor k > 5, the population dynamics converges up to some \c{ct) and, moreover, we clearly 
observe the square root singularity in (n) referenced to this Ac. This suggests that we are detecting 
a bulk transition in the lattice gas, from which we determine ccc by solving Xc{c() = 1/2^, see Eig.[^ 


Appendix D: Improved bounds using the critical exponents 

We can rewrite R(kei H) > exp(n/) = where Pe// = exp(/). In the large n limit, there 
exists a (non-analytic) power series expansion of /(A) near Ac. We learn a few interesting facts 
from the existence of this power series expansion (and in particular, that / does not diverge at the 
critical point). Eirst, Peff 7^ 0 at the critical point. Second, / is known to have non-analytic terms, 
that depend only on the lattice dimension. Eor example, for all 1-D lattices (chains, ladders, etc.), 
the expansion close to the transition point is /(A) = ao -|- ai/2(-^ ~ Ac)^^^ -|- 0{X — Ac), and for 2-D 
lattices, /(A) = 60 + — Ac)^^*’ -I- 0{X — Ac). 

To conclude, near the critical point, the leading power of the non-analytic part of / depend 
only on the dimension of the lattice involved, and not on the specific details of the lattice. As a 
result, the dimension of the satisfying subspace grows faster than one would expect for 1-D and 
2-D lattices. 




